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Statistical analysis of surface distortions of the 70 meter NASA/JPL antenna , located 
at Goldstone , was performed. The purpose of this analysis is to verify whether deviations 
due to gravity loading can be treated as quasi-random variables with normal distribution. 
Histograms of the RF pathlength error distribution for several antenna elevation positions 
have been generated. The results indicated that deviations from the ideal antenna surface 
are not normally distributed . The observed density distribution for all antenna elevation 
angles is taller and narrower than the normal density , which results in large positive values 
of kurtosis and a significant amount of skewness. The skewness of the distribution changes 
from positive to negative as the antenna elevation changes from zenith to horizon. 


I. Introduction 

The surface of the main reflector of a Cassegrain antenna 
is supposed to be a pertect paraboloid, in practice, however, 
there are systematic and random causes of surface irregulari- 
ties. Ripples on the antenna surface as shown in Fig. 1 are 
basically a product of all gravity, wind and thermal loading 
effects. Additionally, some distortions are introduced at the 
factory during panel manufacturing, and some are caused by 
imperfect panel alignment and setting. 

The gain loss, AG, due to reflector surface imperfections 
can be computed in decibels from: 

AG = 10 log 10 ?? (1) 

where the efficiency tj is approximated by Ruze’s (Ref, 1) 
equation as: 



where the root mean square (rms) of surface deviations is 
defined as one -half the change in the RF pathlength and X is 
the operating wavelength (Refs. 2, 3). Fig. 2 illustrates the 
geometrical relationship, between an arbitrary surface displace- 
ment d , normal deviation rz, and the corresponding RF path- 
length deviation 6 . 

Ruze’s Eq. (2) is derived under the assumption that the sur- 
face deviations 5^ at any surface point i are random variables 
having a Gaussian distribution with a zero mean and a standard 
deviation equal to the rms of the surface deviations of the 
reflector (6 t e7V( 0, rms)). Eqs. (1) and (2) are commonly used 
to estimate antenna gain loss even when surface distortions are 
due to deterministic causes such as gravity. 
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The purpose of this article is to verify whether deviations 
due to gravity can be treated as quasi-random variable with 
normal distribution. In order to eliminate truly random devia- 
tions (due to panel misalignment, wind turbulence and local 
thermal effects as well as measurement uncertainty) surface 
deviations simulated by a mathematical finite element model 
of antenna structure are used. 

II. Data Description 


bilities in Ax = 0.1 millimeter segments) are obtained by divi- 
sion of each class frequency by the total number of points 
TV = 7 1 5 . The mean, m , and standard deviation, a, for the 
particular elevation position have been calculated. Histograms 
with the normal counterparts N(m,o) are shown on Figs. 3, 
4, 5, 6 and 7 for the elevation angles 0, 30, 60, 75, and 90 
degrees respectively. In all five cases, histograms do not follow 
the normal distribution. The departure from normality is veri- 
fied by computing the value of chi square x 2 statistics accord- 
ing to the expression: 


The pathlength errors , due to gravity loading, measured 
from the least square best fit paraboloid were obtained for the 
70 meter antenna from a finite element model generated with 
the JPL-IDEAS structural program. Five different antenna 
elevation positions, 0, 30, 60, 75 and 90 degrees, have been 
analyzed. Rigging position, that is 45 degrees, is excluded from 
analysis, because rms is zero at that elevation. Due to antenna 
symmetry, only one half of the dish, represented by N = 715 
points, was modeled. The rms (root mean square) is calculated 
as: 


rms 
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where N is the number of surface points and surface weight- 
ing factors Wj are proportional to the point tributary areas. 
The surface weighting factors w. are normalized to meet the 
condition: 

T, w i = 1 (4) 

/=i 

Statistical analysis is performed on a new rescaled variable: 

x t = (5) 


where 

a. = \J w . N 


k ( p.-p . ) 2 
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( 6 ) 


where K is the number of classes, p { is the probability of 
occurence of observed data in the /- th interval (x { ,x i + Ax) 
and P. is the probability obtained from a hypothetical normal 
distribution N(m y o): 


x.+Ax 

P. = — — f exp[-(x - m) 2 /2a 2 ] dx (7) 

1 a 


The hypothesis that the observed distribution is normal 
can be rejected with probability higher than 0.99 due to the 
large values of y}/N (see Table 1). The same results can be 
obtained by using GFIT program from IMSL (International 
Mathematical and Statistical Library) which computes the 
goodness of fit by x 2 statistics, but it uses equiprobable 
categories (Ref. 4). 

Since it is possible that the estimate of m and o are incor- 
rect it was checked whether values of x 2 statistics can be 
decreased for other normal distributions with different m and 
a. Actually x 2 was minimized with respect to m and a. This 
resulted in finding a normal distribution which best fits the 
observed data. The results, after the best fit, are presented 
in Table 2. Again it can be concluded that with probability 
higher than 0.98, hypothesis of normality, must be rejected. 


Such a scaling allows direct comparison of the standard 
deviation of X and rms. The mean of X is approximately 
equal to zero (in all cases it was less than 0.013 mm) and 
the standard deviation agrees with the rms with better than 
1% accuracy. 

III. Test of Normality 

For each antenna elevation position, a histogram of relative 
frequency distribution has been generated ( K = 30 classes in 
-1.5 to 1.5 millimeter interval). Relative frequencies (proba- 


IV. Non-Normality Characterization 

Since it was verified that gravity induced path length error 
is not normal for any antenna elevation position it is desirable 
to characterize the observed type of non-normality (Ref. 5). 
The two, skewness and kurtosis, parameters characterize the 
type of non-normality. Using higher moments about the mean: 
(with k> 1) 

m k = wT > (x <- m)k ’ * =2 > 3 ’ 4 < 8 > 
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The coefficient of skewness is estimated according to the 
formula: 

m 3 

g ' ~ (m 2 ) 3/2 (9) 

where m 2 and m 3 are the second and third moments, respec- 
tively. If the sample comes from the normal population g } is 
approximately normally distributed with mean zero and 
standard deviation equal \f~6/N. The skewness measures 
the amount of symmetry around the mean. If low values of 
X are bunched close to the mean but high values extend far 
above the mean, this measure will be positive, since the large 
positive contributions of (X -m) 3 , when X exceeds m, 
will predominate over the smaller negative contributions of 
(. X - m) 3 obtained when X is less than m. In the population 
with negative skewness, the lower tail is the extended one, 
shifted to the left. In the population with positive skewness, 
the upper tail is the extended one, shifted to the right. The 
results of skewness for all elevation positions are shown in 
Table 3. The largest negative skewness is confirmed for the 
antenna at the horizon. Skewness decreases when the antenna 
approaches the rigging position (45 degrees). The largest posi- 
tive skewness occurs for 60 degrees elevation and it decreases 
when the antenna approaches zenith position. The same trend 
can be observed on the histograms of Figs. 3, 4, 5, 6, and 7. 

The degree of flattening of a probability density function 
near its mean is measured by the coefficient of kurtosis com- 
puted according to: 



In large samples, g 2 is normally distributed with mean zero 
and standard deviation equal V 24/N. Positive values indicate 
that the probability density function is more sharply peaked, 
that is taller and narrower, than the normal density. Flat- 
topped distributions show negative kurtosis. The kurtosis 
parameter for all antenna elevation positions is presented in 
Table 3. Since the sample’s values of g 2 are much larger than 
the standard error \/24/N , the positive kurtosis parameter is 
confirmed for every antenna elevation position. Again, this 
result of exceeding the normal peak is noticeable on the 
histograms of Figs. 3, 4, 5, 6, and 7. 


V. Summary 

The statistical analysis of the 70 meter antenna surface dis- 
tortions indicated that deviations from ideal paraboloid are 
not normally distributed. The distribution for all elevation 
angles is taller and narrower than normal density, which results 
in large positive values of the kurtosis and a significant amount 
of skewness. The skewness of the distribution changes from 
negative to positive as the antenna elevation angle changes 
from horizon to zenith. Therefore, the gravity induced defor- 
mations violate the assumption of normality in Ruze’s equa- 
tion. It is reasonable, however, to assume that for small depar- 
tures from normality Ruze’s equation will provide a good esti- 
mate of antenna gain loss. Nevertheless, the following open 
problem is posed: what modifications, if any, should be 
applied to Ruze’s equation to incorporate non-zero skewness 
and kurtosis as a function of antenna elevation positions? 
Therefore Ruze’s formula should be verified for non-normal 
quasi-random distribution by comparing it with results pro- 
duced by the exact expression for antenna efficiency in the 
theory of electromagnetic scattering and diffraction (Ref. 6). 
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Table 1. Chi square goodness of fit test: m and o- one obtained 
with standard estimators 


Elevation 

angle 

m 

o 

x 2 / N 

0 

-0.013 

0.646 

0.188 

30 

-0.003 

0.204 

0.600 

60 

0.002 

0.190 

0.239 

75 

0.002 

0.367 

0.183 

90 

0.001 

0.534 

0.090 


Table 2. Chi square goodness of fit test: m and <t minimize x 2 


Elevation 

angle 

m 

a 

x 2 IN 

0 

0.059 

0.576 

0.145 

30 

-0.013 

0.239 

0.233 

60 

0.008 

0.215 

0.156 

75 

-0.010 

0.407 

0.134 

90 

-0.003 

0.474 

0.069 


Table 3. Skewness and kurtosis test 


Elevation 

angle 

Skewness 

Kurtosis 

0 

-1.338 

2.77 

30 

-0.972 

1.82 

60 

0.568 

1.61 

75 

0.429 

2.16 

90 

0.363 

2.99 


0.092 


*JWn 


0.183 
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ig. 5. Histogram of RF path length error distribution at 60 degree 
eievation 
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Fig. 7. Histogram of RF pathlength error distribution at 90 degree 

eievation 




Fig. 6. Histogram of RF pathiength error distribution at 75 degree 

elevation 
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